home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 005 / miscstat.arc / LINREG.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1985-10-07  |  1.0 KB  |  46 lines
  1. 1  'LINear REGression
  2. 2  KEY OFF:SCREEN 0,1,0:WIDTH 40:COLOR 15,1,9:CLS
  3. 10   PRINT "LINEAR REGRESSION"
  4. 20   PRINT 
  5. 30   PRINT "NUMBER OF KNOWN POINTS";
  6. 40   INPUT N:IF N=0 THEN 390
  7. 50  J = 0
  8. 60  K = 0
  9. 70  L = 0
  10. 80  M = 0
  11. 90  R2 = 0
  12. 100   FOR I = 1 TO N
  13. 110   PRINT "X,Y OF POINT";I;
  14. 120   INPUT X,Y
  15. 130  J = J + X
  16. 140  K = K + Y
  17. 150  L = L + X ^ 2
  18. 160  M = M + Y ^ 2
  19. 170  R2 = R2 + X * Y
  20. 180   NEXT I
  21. 190  B = (N * R2 - K * J) / (N * L - J ^ 2)
  22. 200  A = (K - B * J) / N
  23. 210   PRINT 
  24. 220   PRINT "F(X) = ";A;"+ (";B;"* X )
  25. 230  J = B * (R2 - J * K / N)
  26. 240  M = M - K ^ 2 / N
  27. 250  K = M - J
  28. 260   PRINT 
  29. 270  R2 = J / M
  30. 280   PRINT "COEFFICIENT OF DETERMINATION (R^2): ";R2
  31. 290   PRINT "COEFFICIENT OF CORRELATION: " SQR (R2)
  32. 300   PRINT "STANDARD ERROR OF ESTIMATE: " SQR ( ABS (K / (N - 2)))
  33. 310   PRINT 
  34. 320   PRINT "INTERPOLATION: (ENTER X=0 TO END PROGRAM)"
  35. 330   INPUT "X = ",X
  36. 350   IF X = 0 THEN 390
  37. 360   PRINT "Y = ";A + B * X
  38. 370   PRINT 
  39. 380   GOTO 330
  40. 390  PRINT:LOCATE 24,1:PRINT"PRESS <SPACE> TO RUN AGAIN";
  41. 400  LOCATE 25,1:PRINT"PRESS <ESC> TO EXIT";
  42. 410  K$=INKEY$:IF K$="" THEN 410
  43. 420  IF K$=" " THEN RUN
  44. 430  IF K$<>CHR$(27) THEN 410
  45. 10000  RUN"MENU
  46.